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Abstract 

Anomalously strong finite-size effects have been observed for the mean square radius of gyration 
R 2 K of Gaussian random polygons with a fixed knot if as a function of the number N of polygonal 
nodes. Through computer simulations with N < 2000, we find for several knots that the gyration 
radius R? K can be approximated by a power law: R? K ~ N 2u i< , where the effective exponents z/|? 
for the knots are larger than 0.5 and less than 0.6. A crossover occurs for the gyration radius of 
the trivial knot, when N is roughly equal to the characteristic length N c of random knotting. For 
the asymptotic behavior of R\, however, we find that it is consistent with the standard one with 
the scaling exponent 0.5. Thus, although the strong finite-size effects of R\ remain effective for 
extremely large values of N, they can be matched with the asymptotic behavior of random walks. 
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I. INTRODUCTION 



Topological effects on statistical and dynamical properties of ring polymers should be 
quite nontrivial. The topological state of a ring polymer is described by a knot type, and it 
is invariant after its synthesis. Knotted ring polymers or knotted DNAs have been discussed 
theoretically since the 1960s [1], [2], and recently they are synthesized in several experiments 
|3|, |J. Various topological effects on ring polymers have been explicitly studied through 
numerical simulations of random polygons under topological constraints [(], [7], §]. However, 
many questions still remain unsolved even on the average size of a knotted ring polymer in 
solution, which should be the most fundamental quantity in the physics of ring polymers. 



Recently, it has been suggested 0, [10], |TT[] that the average size Rtnv of a random polygon 
with the trivial knot should scale as N u with respect to the number N of polygonal nodes, 
where the exponent v is given by the exponent z/saw of the asymptotic scaling-behavior 
of the the self-avoiding walk (SAW) where z/saw ~ 0.588. Here we remark that random 
polygons correspond to ring polymers with no excluded volume. There has also been a 
conjecture that the effect of topology on the average size of ring polymers could play a 
similar role as the excluded- volume effect, since the topological constraint should effectively 
lead to an entropic repulsion among the monomer segments. Here, we note that the trivial 
knot (or the unknot) is such a knot that is equivalent to an unknotted circle. 

The conjecture on the entropic repulsion seems to be quite interesting, and the anomalous 
scaling behavior with an enhanced exponent should be effective at least for some numeri- 
cal simulations. However, it is not trivial to understand the consequence that the scaling 
exponent should be enhanced and given by that of SAW. First of all, under no topological 
constraint, the average size R of a random polygon scales as N UWN where z/rw = 0.5. Further- 
more, it is not clear whether the scaling exponents of nontrivial knots should be enhanced 
similarly as that of the trivial knot. If they might have the same exponent of SAW, then 
why the total average over all knot types can have the exponent z/rw of random walk ? The 
purpose of this paper is to discuss these questions explicitly through numerical simulations. 
We evaluate the mean square radius of gyration R\ of Gaussian random polygons with a 
fixed knot type K. Discussing the iV-dependence of R\ for several different knot types, we 
show that the anomalous scaling behavior should be considered as a strong finite-size effect 
which could be valid for very large values of N such as 2000. 



2 



We now review some relevant results on the topological effects of ring polymers. Let 
us take a model of random polygons of iV nodes || [7J, which describes ring polymers 
consisting of N Kuhn units at the theta condition. We denote by Pk{N) the probability 
of a given configuration of the random polygon of iV nodes having a fixed knot type K. 
For the trivial knot, it was numerically shown |7|, [5], |T2| that the probability is given by an 
exponential function of N: P tr i V {N) = exp(—N/N c ). For nontrivial knots, the probability 
is well described by the following function of N: P K (N) = C K (N/N c ) m( - K ^ exp(-N/N c ), 
where we call N c and m(K) the characteristic length of random knotting and the topological 
exponent of the knot, respectively ||13|| . The value of N c is model-dependent, and is roughly 
given by 340 for the Gaussian random polygon || 13]. We remark that the number N c is 



important in the analysis of topological effects with the blob picture [11 



The mean-squared gyration radius R 2 K under the topological constraint of knot K has 



been discussed for some models of self-avoiding polygons (SAPs) in Refs. |14], [15], [L6|, [17], [18|, 
19| , |2"0fl . In the lattice model, it is shown that the asymptotic behavior of R 2 K is consistent 



with that of the RG theory where in the large N limit the ratio R K /R comes close to 1.0 



for any knot. However, for the cylinder model of SAPs [21|, it is found [Kj that the limit of 
the ratio depends on the cylinder radius which controls the excluded-volume. For a lattice 
model of random polygons [^] , R\ has been evaluated for the trivial and trefoil knots with 
small polygons of N < 200. 



II. METHODS OF SIMULATIONS 

We now introduce the Gaussian random polygon |J. Let A 1; • ■ ■ , denote the position 
vectors of the nodes of a configuration of the Gaussian random polygon of N segments, 
and ,un the jump vectors such that Uj = Xj — Xj-i for j = 1, ... ,JV. Then, 

the Gaussian random polygon has the following distribution function of the jump vectors: 
P(ui, ■ ■ ■ ,un) = Const, x exp(— {u\ + • ■ ■ + u%)/2) 5{ui + ■ ■ • + un) ■ Here 5(x) denote 
Dirac's delta function in three dimensions. We construct M configurations of the Gaussian 
random polygon with N nodes by the conditional probability distribution 

P(u f , u u -.., Vi) = (2vr)- 3/2 exp (- ^J^V i + WZJ^i A ^ 
Let us consider M samples of randomly constructed polygons with N nodes. We denote 
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by Mk the number of polygons with knot K. We determine the number Mk by enumerating 
such polygons in the set of M polygons that have the same set of values of the following two 
invariants for knot K: the determinant of knot A#(— 1) and the Vassiliev invariant V2{K) 
of the second degree P31 E3). 



III. ANOMALOUS FINITE-SIZE BEHAVIORS OF R\ 

Let us discuss the numerical data of our simulations. For Gaussian random polygons 
with several different numbers iV up to 1900, numerical estimates of R 2 and R 2 K have been 
obtained for the four knots: the trivial, trefoil (3i) and figure-eight (4i) knots, and the 
composite knot consisting of two trefoil knots (3i(j3i). Here, we take M = 10 5 in all the 
simulations. 

We recall that under no topological constraint, the mean square radius of gyration R 2 of 
a random polygon with N nodes is defined by R 2 = J2nm=i((^n — Rm) 2 ) /2N 2 . Here R n is 
the position vector of the n-th node and the symbol (■) denotes the statistical average, which 
is given by the average over M polygons in the simulations. For a knot K, the quantity 
R 2 K is given by R\ = J2i=a RkJMk, where R 2 Ki denotes the gyration radius of the z-th 
Gaussian random polygon that has the knot type K, in the set of M polygons. In terms of 
R 2 K , R 2 is given by R 2 = J2 K M K R 2 K /M. 

In Fig. 1, the graphs of the ratio R 2 K jR 2 against the the number N are depicted for the 
trivial, trefoil and figure-eight knots. The fitting curves in Fig. 1 are given by R 2 K /R 2 = 
TkN 2AuC k . The best estimates of the fitting parameters are given in Table 1. 

We see in Fig. 1 that the ratio R\jR 2 is not constant with respect to N. The ratio 
R\jR 2 increases monotonically for all the knots. For the trivial knot, the ratio R 2 riv /R 2 is 
always larger than 1.0. On the other hand, for the other three knots (K=3\, 4i, 3ijj3i), the 
ratio R 2 K jR 2 is smaller than 1.0 when N is small. However, the ratio R 2 K jR 2 becomes larger 
than 1.0 when N is large. Thus, the topological constraint gives an effective swelling for the 
case of large N. Here we have a conjecture that for any nontrivial knot K, R 2 K should be 
larger than R 2 if N is large enough. 

Let us consider the plot of the trivial knot shown in Fig. 1. There is a nontrivial finite-size 
behavior: the ratio R 2 riv /R 2 is constant wit respect to N when N is very small, while when 
iV > N c , it can be approximated by a scaling behavior as R 2 riv /R 2 ~ N 2u ^, at least up to 
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N = 2000. This crossover phenomenon should be consistent with the recent theory given 
by Grosberg JO]. However, the effective exponent vff iv is much smaller than the exponent 
^saw with respect to the errors. In fact, we have the numerical estimate: v^ iv ~ 0.545. 

For the case of nontrivial knots (3i, 4 1; 3i(j3i), the ratio R 2 K / R 2 is well approximated by 
the power law: R 2 K /R 2 « T k N 2Au k for the range from N = 100 to iV = 2000. Furthermore, 
there is no crossover for the nontrivial knots: we do not find any change in the slope of the 
graphs near iV ~ N c . It is also remarkable from Table 1 that the exponent Az/|? strongly 
depends on the knot type. In particular, the effective scaling exponent of the composite 
knot 3i#3i is almost as large as the exponent z^saw, while that of the trefoil knot is given 
by 0.561, which is rather smaller than z/gAw with respect to the errors. 

Let us consider the three fitting lines of Fig. 1. Then, we see that the three lines 
become very close to each other at around N = 2000. Furthermore, it is suggested from 
the simulations that when N becomes close to 2000, the values of R 2 K for the four knots 
(if=trivial, 3i, 4i and 3i(j3i) should become almost equal to each other. In fact, up to 
N = 1900, the values of R 2 K for the nontrivial knots are always smaller than or equal to 
that of R 2 riv in our simulations. If the power-law approximation might be valid also for 
iV > 2000, then R 2 riv would become much smaller than R 2 K of the three nontrivial knots for 
large N and it would be inconsistent with the numerical results obtained so far. Thus, we 
may conclude that the approximation of R 2 K with the power law should be valid only when 
N < 2000. Therefore, in order to study the iV-dependence of R 2 K for iV > 2000, we need 
another independent analysis. 



IV. ASYMPTOTIC BEHAVIOR OF R 2 K 

Let us discuss the asymptotic behavior of R 2 K . When N is very large, we may assume 
the following expansion: R\ = A k N 2uk (1 + B K N- A + 0(1/ N)). It is consistent with 
renormalization group arguments, and hence it should be valid when N is asymptotically 
large. For the ratio R 2 K /R 2 we have the following expansion: 

R\jR 2 = (A K /A)N 2Au *{l + (B K - B)N~ A + 0(1/ N)) , (2) 

where Auk = vk — ^rw- Here we recall that z/rw — 0.5. For each of the four knots, we 
have applied the formula (||D to the 10 data points with iV > 1000 shown in Fig. 1. Here, 
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when we assign the condition of N > 1000, we have taken into account the strong finite-size 
effects of R 2 K such as the crossover of the trivial knot. Another recent study [^5j on the 



cylinder model of SAPs |20| , [2l|j shows that the asymptotic scaling behavior is seen only 
when N > 1000 for the case of very thin cylinders with the cylinder radius r = 0.001. 

The best estimates of the fitting parameters of the formula (0) and the x 2 values are 
listed in Table 2. From the results, we may conclude that the asymptotic expansion is 
consistent with the numerical values of R 2 K for iV > 1000. It is remarked that the \ 2 values 
in Table 2 are less than 10 for the four knots. Moreover, the best estimates are compatible 
with several different viewpoints. For instance, the estimate of IAvk is given by about 0.03 
and independent of the knot type. This leads to an estimate of the exponent: vk rs 0.515. 
The value can be considered as equivalent to the exponent ^rw ; with respect to the errors 
of the analysis. The fact that the exponent vk is independent of the knot type is consistent 
with the interpretation on the lattice model of Refs. |17| , [Tj|. The estimated values of the 
amplitude ratio Ak/A for the four knots also seem to be independent of the knot type. We 
note that we have assumed A=0.5 in constructing Table 2. When we set A=1.0, similar 
values are obtained for the best estimates of the fitting parameters. 

Let us consider a formula which effectively describes the iV-dependence of R 2 K for iV > 
2000. Assuming Avk — in the asymptotic formula (^), we have the following: R 2 K /R 2 = 
ozk(1 + Pk N~ a + 0(1/N)). Here we have replaced by olk and (3k, A k /A and B K — B, 
respectively. Applying the formula to the numerical data of R 2 K of the Gaussian random 
polygon for iV > 1000, we see that it gives good fitting curves to the data. The best 
estimates of the parameters are shown in Table 3. Interestingly, they are rather close to 
the best estimates for the cylinder model of SAPs with a very small cylinder radius, which 
are obtained by applying the same formula to the data of R 2 K in Ref. [pOR . In Table 3, the 
parameter olk is roughly given by 1.5 for the Gaussian random polygon. On the other hand, 
we have the similar value for the cylinder model with the cylinder radius r = 0.001 as shown 
in Fig. 3 of Ref. p0|| . We also find in Table 3 that olk ~ 1-5 for the four knots. It follows 
that the mean size Rk of random polygons with a specified knot K, such as the trivial, 3i, 
4i, and 3i#3i knots, is larger than the average size R of random polygons over all knots 
in the asymptotic regime. However, it is consistent with the observation in Fig. 1 that the 
ratio R 2 K /R 2 increases monotonically and approaches 1.3 or 1.4 when N ~ 2000 (see also 
Ref. [EDI). 
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V. CONCLUSION 



We have shown that R 2 K of Gaussian random polygons have strong finite-size effects which 
should be valid for extremely large values of N such as N = 2000. If we remove the finite-size 
effects from the data analysis, then the asymptotic behavior of R\ is given by the standard 
one with the critical exponent u-^w of random walks. Here we note that the main result 
should be valid also for the numerical data of other models of random polygons |26| . 
Thus, the studies || [TU], [TT| associated with the conjecture of the effective entropic-repulsion 
should be important for describing the strong finite-size effects of random polygons, which 
could appear practically in any system of ring polymers in solution at the theta condition. 
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TABLE I: Best estimates for the fitting lines in Fig.l describing the anomalous scaling behavior: 
Rl/B 2 = T K N 2Av f. Here Avf = - z/ RW . For the trivial knot, the fit is obtained from the 
data points with N > 400. For the trefoil (3i) and figure-eight (4i) knots, and the composite knot 
of 3i#3i, the fitting lines are obtained from the data for N > 100. 



knot type T K 2Avf \ 2 



triv/ave 


0.600±0.012 


0.090±0.003 


6 


tre/ave 


0.514±0.004 


0.121±0.001 


14 


4i / ave 


0.431±0.006 


0.143±0.002 


11 


3i(J3i/ave 


0.398±0.005 


0.153±0.002 


25 


TABLE II: 


Best estimates of the fitting parameters of the asymptotic formula: R 2 K /R 2 = 


{A K /A)N 2A 


VK (l + (B K - B)N- A ). Here Au K = v K - z/ RW . 




knot type 


A K /A 


B K -B 2Av K 


x 2 


triv/ave 


1.163±3.22 


-3.746±17.609 0.027±0.309 


3 


tre/ave 


1.131±1.837 


-5.333±9.486 0.033±0.183 


9 


4i / ave 


1.129±3.837 


-5.860±19.170 0.033±0.385 


7 


3ijj3i/ave 


1.166±1.662 


-5.928±8.071 0.028±0.161 


9 


TABLE III: Best estimates of the fittinj 


I parameters of the formula effectively describing the large- 


N behavior of the ratio: R 2 K /R 2 = ax 


(1 + P K N~ A ). 




knot type 


a K 


Pk 


x 2 


triv/ave 


1.472±0.042 


-5.187±0.833 


3 


tre/ave 


1.516±0.026 


-6.923±0.483 


10 


4i / ave 


1.507±0.0056 


-7.375±1.1010 


7 


3it)3i/ave 


1.493±0.024 


-7.270±0.446 


9 
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FIG. 1: Logarithmic plot of the ratio R 2 K /R 2 versus the number N of polygonal nodes of the 
Gaussian random polygon for the range from N = 100 to N = 1900. Numerical estimates of 
R 2 K /R 2 for the trivial, trefoil (31) and figure-eight (41) knots are shown by black circles, black 
squares and black triangles, respectively. In the inset, the enlarged figure shows the logarithmic 
plot of the estimates of R 2 K /R 2 from N = 1000 to 1900 for the trivial, trefoil and figure-eight knots. 
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